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Self-organization of a long-lived structure is one of the remarkable characteristics of macroscopic 
systems governed by long-range interactions. In a homogeneous magnetic field, a non-neutral plasma 
creates a “thermal equilibrium” which is a Boltzmann distribution on a rigidly rotating frame. Here, 
we study how a non-neutral plasma self-organizes in inhomogeneous magnetic field; as a typical 
system we consider a dipole magnetic field. In this generalized setting, the plasma exhibits its fun¬ 
damental mechanism that determines the relaxed state. The scale hierarchy of adiabatic invariants 
is the determinant; the Boltzmann distribution under the topological constraint by the robust adi¬ 
abatic invariants (hence, the homogeneous distribution with respect to the fragile invariant) is the 
relevant relaxed state, which turns out to be a rigidly rotating clump of particles (just same as in a 
homogeneous magnetic held), while the density is no longer homogeneous. 
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INTRODUCTION 

Self-organization of a long-lived structure (inhomo- 
geneity of physical quantities) is often observed in macro¬ 
scopic systems governed by long-range interactions such 
as gravity (creating astronomical systems like galax¬ 
ies [T]), electromagnetic force (creating plasma systems 
like magnetospheres im or particle traps m , or mag¬ 
netic interaction (Hamiltonian mean-held systems mod¬ 
eling magnetism [6l[7])- The common physics is described 
by the Vlasov equation coupled with a relevant held equa¬ 
tion. The long-lived structure is a particular stationary 
solution that is “robust” against microscopic perturba¬ 
tions. Here, we put the problem into the perspective 
of non-canonical Hamiltonian mechanics [8] , and show 
that the self-organization occurs on a leaf of topologi¬ 
cally constrained phase space; the topological constraint 
originates from the adiabatic invariants, which dehnes a 
macroscopic hierarchy [9] . As a specihc system, we con¬ 
sider a non-neutral (single species) plasma in a dipole 
magnetic held. Let us begin by explaining how this prob¬ 
lem is interesting from both basic and applied physics 
viewpoints. 

When a non-neutral plasma is put in a homogeneous 
longitudinal magnetic held, it is spontaneously conhned, 
relaxing into a “thermal equilibrium” on a rigidly ro¬ 
tating frame laiioiiii]. Canceling the self electric held 
by the Lorentz-transformed electric hied, the Boltzmann 
distribution yields a homogeneous density prohle inside 
the conhnement region. The constant density prohle and 
the constant angular momentum prohle are the simulta¬ 
neous characteristics of the relaxed state] there remains 
no free energy to excite macroscopic perturbations (as far 
as the system conserves the total angular momentum and 
is isolated from other energy sources like other species of 
particles or external electromagnetic helds). However, 
the consistency of the density distribution (dictated by 
the statistical mechanics of particles) and the electric po¬ 


tential (dictated by the held equation) relies heavily on 
the specialty of the homogeneous longitudinal magnetic 
held. Here, we investigate whether such a relaxed state 
exists in an inhomogeneous magnetic held. Experimen¬ 
tally it does exist in a dipole magnetic held [mils]; the 
charged particles self-organize a rigidly rotating clump. 
However, the density is no longer homogeneous. The 
aim of this study is to reveal the underlying principle 
that governs generalized relaxed states. Conhnement of 
charged particles (especially antimatter particles) in a 
toroidal magnetic bottle has many advantages, for ex¬ 
ample, making possible to conhne high-energy particles 
produced by isotopes or accelerators, or to conhne differ¬ 
ent spices of positive and negative charges simultaneously 

mHin]. 

Here we invoke the theory of phase-space foliation (or, 
topological constraint), and dehne a relaxed state as 
a thermal equilibrium on a leaf of phase space [9] . In 
the present argument, the adiabatic invariants of magne¬ 
tized particles embody such foliated phase space. In an 
axisymmetric magnetic held, magnetized particles have 
three different adiabatic invariants, i.e., the magnetic mo¬ 
ment /i, the action Jy of bounce motion, and the ac¬ 
tion (canonical angular momentum) Pq of the toroidal 
drift [17]; see Fig.[^ We may approximate Pe/q {q is 
the charge) by the magnetic hux function ip such that 
B = X \/0 {0 is the toroidal angle). When the mag¬ 
netic held is sufhciently strong, the corresponding fre¬ 
quencies dehne a hierarchy: uJc (cyclotron frequency) ^ 
ujb (bounce frequency) 7^ Ud (drift frequency). Hence, 
is the most fragile constant —the homogenization with 
respect to ip yields the relaxed state on the hrst (macro¬ 
scopic) hierarchy of the adiabatic invariants. Needless to 
say, the ultimate relaxed state is achieved at the maxi¬ 
mum homogeneity (maximum entropy) after destroying 
all adiabatic invariants, and it is the thermal death. 

We dehne the “relaxed state” by a distribution func- 
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initial 



The action Jy is conjugate to 13^: 

. _ dH,, _ ^11 dH,, 
'' IT dC 


For the periodic bounce motion, f{dHgc/d'3b)d'3}y = 
f dHgc = 0. Integrating ^ over the cycle of bounce mo¬ 
tion yields the bounce-average (Jy) = constant. When 
we calculate macroscopic quantities (like the total en¬ 
ergy or the total action), we evaluate Jy as the adiabatic 
invariant (Ty). 

The drift frequency (including all grad-B, curva¬ 
ture, and ExB drifts) is given by bounce-averaging the 
toroidal angular velocity 


ujd = 0 


dHgc _ d^ ^ 

d^p ^ d^ dtjj ^ dt^ 


(4) 


FIG. 1: A typical orbit of a magnetized particle in a dipole mag¬ 
netic field, which has a hierarchy of three different frequencies. Here 
we assume the parameters of the RT-1 device CU; the magnetic 
field at the confinement region is about 0.05T. For a particle having 
an isotropic energy of 50eV, ouc ~ IGHz, ~ IMHz, and ~ 
IkHz. 

tion / that has no -0 dependence, i.e.. 


For the equilibrium to be non-trivial, we demand that 
the total canonical angular momentum f ipfd^z to be a 
non-zero constant {d^z denotes the volume element of the 
phase space). 


In a homogeneous magnetic field, both ujc and are 
constant, and then Q evaluates the ExB drift frequency. 

In terms of the constants of motion Hgc, Jy, and -0, a 
general equilibrium solution of the drift kinetic equation 
(such that {Hgc^f} = 0 ) is written as /(i^^c, J||, V^)- 

The relaxed state is the special solution that maximizes 
the entropy S = — f f log f d^z under the constraints on 

1. the total particle number N = f f d^z^ 

2 . the total energy E = J Hgcf d^z^ 

3. the total magnetic moment = f (if d^z^ 

4. the total bounce action = f J\\f d^z^ 

5. the total angular momentum = J ^pf d^z. 


KINETIC MODEL OF MACROSCOPIC 
RELAXED SATE 


The variational principle yields 




In order to formulate the model with taking into ac¬ 
count the hierarchy of adiabatic invariants, we write the 
Hamiltonian of a particle as 

Hgc = LdclJi + uJhJ\\ + #• ( 2 ) 

We have omitted the kinetic energy {Po — qpj)‘^/ (2mr^) of 
the toroidal drift velocity by approximating Pq = qp; [18] . 
By the symmetry, the toroidal angle 0 is not included in 
Hgc. The gyro angle ^c (which is conjugate to the mag¬ 
netic moment fi) is coarse-grained by replacing dc with 
cjc, and is completely eliminated from Hgc (i.e., Hgc dic¬ 
tates the motion of the guiding center of the gyrating 
particle). However, the bounce angle ('d^) is not ignored, 
because the frequencies ujc and as well as the elec¬ 
tric potential (j) are functions of the spacial coordinates 
includingHere we choose pj and ( (the parallel co¬ 
ordinate along each magnetic surface, the level-set of pj) 
as the spatial coordinates (then, = 7 rC/^y; ^y is the 
bounce orbit length). 


where Z (normalization factor), P (inverse temperature), 
7 i, 72 , and 73 are constants related to the Lagrange mul¬ 
tipliers on A, E, Cj||, and C- 0 . 

While we derived for a given set of constants, we 
may, alternatively, regard /t as a Boltzmann distribu¬ 
tion on a grand-canonical ensemble parametrized by the 
aforementioned macroscopic quantities, and then, we in¬ 
terpret 7 i, 72 and 73 as the chemical potentials pertinent 
to the changes in the action variables /i, Jy, and -0, re¬ 
spectively (remember the parallel relations between the 
“energy level” and the frequency, as well as between the 
“particle number” and the action variable, in analogy 
with the Landau levels in quantum theory). 

Finally, the determining equation 0 must be satisfied, 
which is equivalent to 

= 73 (= constant). (6) 

Remembering Q, we find that 0 implies rigid rotation. 
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The relaxed-state plasma occupies a finite domain that 
is surrounded by a magnetic surface, i.e., we can connect 
the distribution function /t to the vacuum / = 0 at 
some level-set 7/^ = t/;*. Invoking Heaviside’s step function 
y(' 0 * — (which is zero inside the plasma region), we 
may write the extended distribution function as 

= lim ■ fr, ( 7 ) 


which is in equilibrium {{Hgc^ fT} = 0 ) and relaxed: 
dfr/di^ = 0 except at the boundary (the boundary is 
not “relaxed” as it reflects the constraint given by 
later, we will show how the boundary is determined). 

Neglecting the current generated by the plasma, is a 
given function. Then, may be viewed as a determin¬ 
ing equation for the electric potential 0 (in addition to 
the term Jy depends on 0 in a complex way; however, 
if the electric field is much stronger than the thermal en¬ 
ergy, we may approximate B • ~ 0 , and then, Jy is 

independent to 0 , and we may put (j) = 0 (' 0 )). 

The electric potential (j) included in Hgc must be con¬ 
sistent to the field equation 

—V^(/) = 47 rgp, ( 8 a) 

!•= j hM<i‘v = j ( 8 b) 


The existence of a self-consistent field (j) satisfying both 
© and ( [ 8 a[ ) is not at all obvious. In what follows, we 
will construct non-trivial solutions; one is the well-known 
“thermal equilibrium” in a straight homogeneous mag¬ 
netic field, and the other is a new solution (numerical) in 
a dipole magnetic field. 


THERMAL EQUILIBRIUM IN A 
HOMOGENEOUS MAGNETIC FIELD 


RELAXED STATE IN A DIPOLE MAGNETIC 
FIELD 


In a dipole magnetic field, both ujc and vary as 
functions of the spacial coordinates and Because of 


the inhomogeneous Jacobian weight in the integral ( 8 b), 


the parameters 71 and 72, included in /t, cause a change 
in the profile of p nng, which is in marked contrast to 
the foregoing case of homogeneous magnetic field. 

Assuming some (j) in /t^, we calculate p by ( 8 b), and 
then solve ( 8 a) for a new (j). Iterating this process, we 


obtain a self-consistent (j) and a kinetic equilibrium /t- 
Next we have to vary the parameters 71, 72, and 73 to 
find the reiaxed state. 

Since the ieft-hand side of ([^ contains p and Jy, we 
cannot satisfy for each partrSe having different p and 
Jy (excepting the case when duoc/di^ = du^yjd^ = 0 , as it 
is in a homogeneous magnetic fieid). Instead, we demand 
that the macroscopic drift veiocity 


ve = p J vefr^^^^dpdvii 


( 10 ) 


is rigid rotation {vq is directiy evaiuated by the orbit cai- 
cuiations). The chemicai potentiais 71,72, and 73 are 
the controi parameters to be optimized to yieid rigid ro¬ 
tation. 

Figurej^ shows an exampie of soiution. Here we as¬ 
sume parameters that simuiate the RT -1 experiment m- 
The temperature /3 ^ = 50 eV is chosen to be the typi¬ 
cal energy of the injected electrons. Other parameters 
are (71,72,73) = (- 8.0 x T,- 7.2 x 10 S-\ 1.4 x 
10 "^ Wb~^ eV), by which /^yi/i, /372Jy, and / 3 g 0 

are all of order ^ 1 (implying that all terms in /t 
play essential roles in characterizing the relaxed state). 


First, we put the classical solution into the new per¬ 
spective formulated here. In a homogeneous longitudinal 
magnetic field B = BqS/z {Bq is a constant), uJc and 00^ 
are constants. We may put 1/^ = Hor^/2 and 0 = 0 (' 0 ) 
(assuming that the relaxed state is homogeneous with 
respect to C = 2:). The relaxed-state condition 0 reads 
qd (\)/ d ^ — 73, which yields = (73/^)'^, and 

fj. = ^-lg-/5[(^c-7i)M+(^b-72)J||]^ 

Since this distribution function has no spatial depen¬ 
dence, p = po (constant). On the field equation ( [8a| ), 
the potential 0 = (73/^)^^ = (73/g')Hor^/2 is consistent 
to the flat density po, if 73 = — 27 rg^po/^o = —^^/(^ccjc). 
The parameter 73 controls the density po (while 71 and 
72 change the velocity anisotropy). Given the particle 
number N (per unit length in z), po = A/(7rR^), where 
R is the radius of the plasma column. The total angular 
momentum evaluates Cp = BoNR‘^/A. Hence, for given 
N and Cp, we obtain R = /BoN. 


NUMERICAL ANALYSIS 


Here we study how general equilibrium solutions (con¬ 
sistent 0 and /t satisfying /t} = 0 and the Poisson 

equation ( [8a| ) simultaneously) vary as the parameters are 
changed, and how these parameters can be optimized to 
send the equilibrium to the relaxed state. 

For the convenience, we introduce indexes to evaluate 
the “relaxation”: 



(lla) 

(llb) 

(llc) 


i.e., ujd the spatially averaged toroidal drift frequency, 
the associated standard deviation, and y a measure 
of rigidity of the rotation. 
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FIG. 2: A relaxed state in a dipole magnetic field: (a) self- 
consistent density profile p, (b) the distribution of the angular 
velocity ojd, and (c) the contours of the electric potential 0 
(black lines) and the magnetic flux function ?/; (the magnetic 
field lines; red lines). 



FIG. 3 : The rigidity parameter y as functions of 71, 72, 73 
and 00- Each parameter is normalized by the correspond¬ 
ing central value (71,72,73,^0) = (—8.0 x 10“^T, — 7.2 x 
10 ® 1.4 X 10 ^ Wb~^ eV ,—250 V), and changed indepen¬ 

dently. 


p). Here, other parameters are fixed at = —250 
= 50eV, and ( 72 , 73 ) = (-0.72 x 10^5-1,1.25 x 
10 ^ Wb~^ eV). As 7 i is increased, the density approaches 
to the distribution p oc H [4]. However, the rigidity of the 
rotation is a weak function of 71 (see Fig. [^. 


Ki = -0.9 10 


oJciikHz) 



0.3 0.5 0.7 0.9 


= -2.4 10-1 


Yi 


-2.9 10-1 7 


In Fig. the behavior of x as a function of the control 
parameters 71 , 72 , 73 , and 0 o (the electric potential at 
the coil surface) is shown. A detailed explanation for 
each parameter is given in the following subsections. 


FIG. 4 : Plots of the angular frequency uJd{kHz) of rotation 
for four different values of 71 (T). 


Chemical potential 71 of the magnetic moment p 

The dependence of the equilibrium state on the pa¬ 
rameter 7 i (the chemical potential of the magnetic mo¬ 
ment p) is shown in Figs. (plots of uud) and (plots of 


Chemical potential 72 of the bounce action Jy 

The chemical potential 72 of the bounce action Jy 
has a rather strong influence on the equilibrium distri¬ 
bution; see Figs. (plots of uud) and (plots of p). 
With the optimum value 72 = 0.7 x 10 ^ 5 “^, the den- 
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FIG. 5 : Plots of the density p (normalized) for four different FIG. 7 : Plots of the density p (normalized) for four different 
values of 71 (T). values of 72(5”^). 


Chemical potential 73 of the angular momentum 


sity p has a broad distribution with a highly homo¬ 
geneous drift frequency However, when 72 is in¬ 

creased, the density profile shrinks to a disk-like shape 
with a strong shear in Here, other parameters are 
fixed at = —250H, = 50eV, and ( 71 , 73 ) = 

(-0.795 X 10-^ T, 1.25 x lO^VF^-^eH). The change of 
X as a function of 72 is given in Fig. 


cOdikHz) 

72 = —0.7 10^5“^ 72 = “4.3 10^5“^ 
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The chemical potential 73 of the angular momentum -0 
also has a strong influence on the equilibrium; see Figs. 
(plots of ujd) and[^ (plots of p). As 73 is increased, the dis¬ 
tribution function becomes / ~ g^j^d then the 

drift frequency scales as uJd ^ since uJd = qd(\)/d^ 
and, at z = 0 , ^ and 0 ^ (if we assume 

a constant particle number per flux tube volume m, 
p r ^ which gives Acj) ^ r ^). Here, other pa¬ 
rameters are fixed at 0o = —250 H, = 50eV, and 
( 71 , 72 ) = (-0.795x10-^ T,-0.72x10^5-1). The change 
of X as a function of 73 can be found in Fig. 


oidikHz') 

K3 = 3.5 10^ Wb-'^eV = 1.85 10^ Wb-^eV 



0.3 0.5 0.7 


73 = 3.35 10^ Wh-^eV 


73 = 4.85 1 0^ Wh-^eV 


FIG. 6 : Plots of the angular frequency uJaikHz) of rotation 
for four different values of 72 (s ^). 


FIG. 8 : Plots of the angular frequency uJd{kHz) of rotation 
for four different values of ^3{Wh~^ eV). 
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(ji)a,(kHz) 

(Po = -200V (Po = -2S0V 



(Po = -3S0V 


(Po = -400V 


FIG. 9 : Plots of the density p (normalized) for four different 
values of eV). 


FIG. 10 : Plots of the angular frequency uJd{kHz) of rotation 
for four different values of 0o(I^)- 


Inner boundary potential 0 o 

In addition to the four parameters /3, 71 , 72 , and 
73 , there is another important parameter 0 o that is the 
boundary value of the potential 0 at the surface of the 
internal magnet. By changing the confinement is 
dramatically improved; see the experimental result re¬ 
ported in [20] . On the levitated magnet experiment [ 12 ] , 
however, is spontaneously determined, because the 
coil surface is floating. Figures and [m respectively, 
show the distribution of the drift frequency and the 
density p for four different values of 0o- Other pa¬ 
rameters are fixed at = 50eV and ( 71 , 72 , 73 ) = 
(-0.795 X 10-1 T,-0.8 x 10 ^ 5 - 1 , 1.35 x lO^VF^-ie^). 
Interestingly, as changes, the confinement region (the 
density clump) moves, while the rigidity of the rotation 
is not infiuenced by 0 o; see FiglU 



FIG. 11 : Plots of the density p (normalized) for four different 
values of 0o(l^)- 


CONCLUSION 

Putting a classical knowledge into a wider perspective, 
a deeper principle may emerge; here we have described 
such an example of paradigm shift in the study of the 
relaxed state (or, thermal equilibrium) of charged parti¬ 
cles. Formulating a relaxed state on a topologically con¬ 
strained macroscopic phase space, we have found that a 
rigidly-rotating equilibrium can self-organize even in an 
inhomogeneous magnetic field. In an experimental sys¬ 
tem, tuning of the parameters 71 , 72 , 73 , and occurs 
spontaneously through the relaxation process. Among 
them, 72 and 73 influence strongly on the profile of 
ujd. Whereas we formulated the equilibrium assuming 


that the corresponding and are given (then, the 
Lagrange multipliers 72 and 73 are determined by pre¬ 
scribed (7j|| and C^), the plasma may change them (by 
dissipating the constants of motion) in order to relax into 
the thermal equilibrium; the self-organization is a process 
that selects optimum and (hence, 72 and 73 ). In 
fact, these two actions are relatively “fragile” with re¬ 
spect to the other constant C^. 

The present model of relaxed states differs from the 
previously formulated Boltzmann distribution of a neu¬ 
tral plasma [4] , which is created by constraining only 
and (7j|| (in addition to the standard constraints N and 
E) in maximizing the entropy; freeing means that 
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is not a conserved quantity. In the present formula¬ 
tion, however, we also constrain while we demand 
df /d^p = 0 diS the criteria of the relaxed state. The lat¬ 
ter is a weaker criterion, i.e., the present solution is not 
necessarily the maximum entropy state of the former set¬ 
ting, as far as is freed from the kinetic energy (i.e., we 
omit the energy of the drift velocity) and is deemed as 
a spatial coordinate variable. This is obvious by putting 
73 (the Lagrange multiplier on = 0 in the thermal 
equilibrium Then, the solution becomes ^ 0, 
p ^ 0 with Pqcp constant (= oc), i.e., the radius of 
the plasma column diverges, implying no confinement. 
The constraint on the total angular momentum yields 
a finite radius of confinement. 
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